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Abstract 

Siver A.S., Ezhela V.V. On the CODATA Recommended Values of the Fundamental Physical Constants: 
V3. 2(1998) & V4.0(2002): IHEP Preprint 2003-34. - Protvino, 2003. - p. 18, tables 12. 

With the help of special program package PAREVAL designed in Mathematica system we reproduce 
values of basic fundamental physical constants obtained by NIST and recommended by CODATA for in- 
ternational usage since 1998. In our adjustment we use the input data and methods published by NIST in 
1998. 

It is shown, that the detected earlier inaccuracy of published by NIST correlations (that made the NIST 
1998&2002 results doubtful) are, most probably, due to inadmissible independent rounding. 

The simple estimate of the critical numbers of decimal digits in the independently rounded correlation 
coefficients is obtained. Further independent rounding of "critically rounded" correlations can lead to the 
non positive semi definite correlation matrices and hence is inadmissable. 

It is demonstrated by a few examples that the poor presentation of the correlated random quantities in 
the scientific literature is a common bad practice and is argued (once again) that the common standard for 
presentation numerical values of correlated quantities in publications and sites is urgently and badly needed. 

AHHOTail,HH 

CnBep A.C., EjKejia B.B. irHCJieHHbix snaHeHHax cjjyHflaMeHTajibHbix (jansH^ecKHx nocTOHHHMx, peKOMeHflOBanHbix 
CODATA BepcHH V3. 2(1998) n V4.0(2002): DpenpHHT Ha)B9 2003-34. - DpoTEHHO, 2003. - 18 c, 12 Ta6ji. 

C noMombK) cneiinajibHoro naxeTa nporpaMM PAREVAL, paspaGoTanHoro b CHCTCMe Mathematica, 
BOcnpoHSBefleHbi SHa^ennfl 6a30Bbix (JjyH^aMeHTajibHbix (J)H3HHecKHx nocTOHHHbix NIST, peKOMeHflOBanHbix 
CODATA fljisi MeiK^ynapo^Horo HcnojibsOBanHfl c 1998 ro^a. HoKasaHO, hto o6Hapy}KeHHbie panee HeTOirHOCTH 
snaHeHHii K034)4)HLi,HeHTOB KoppejiHu,HH norpemHOCTefi KOHCTanT 1998 ro^a o6ycjiOBJieHbi, CKopee Bcero, 
Hx neflonycTHMbiM nesaBHCHMbiM OKpyrjieHneM. Ko3(|)4)HLi;HeHTbi Koppejismnn norpemHOCTeii KOHCTanT 
2002 rofla TaK»ce HcnopHCHbi neflonycTHMbiM OKpyrjieHHCM, h hx nejibSfl ncnojibsOBaTb b BbicoKOTO^Hbix 

Bbl^HCJIBHHflX. 

nojiy^ena npocTaa ou,eHKa KpnTH^ecKoro HHCJia flecHTH^Hbix snaKOB b snaHeHHAx K03(|)(|)Hu,HeHTOB 
KoppejiHu,HH, ^ajibHeiimee nesaBHCHMoe OKpyrjieHne KOTopbix MOJKeT npnBO^HTb k noTepe nojiojKHTejibHoii 
nojiyonpeflejieHHOCTH MaTpHu;bi KoppejiaiiHii. 

Ha HGCKOJibKHx npHMepax noKasano, hto neKoppeKTHoe npeflCTaBjieHne SHa^reHHii KoppejinpoBaHHbix 
cjiyHaHHbix BejiH^HH B ny6jiHKai];Hflx — sto pacnpocTpaneHHaa Bpe^Hafl npaxTHKa. /JeKjiapHpyeTCH Heo6xoflHMOCTb 
pa3pa6oTKH H HcnojibsOBaHHH CTan^apTa npe^CTaBJienHH OLi;eHOK KoppejinpoBanHbix cjiy^aiiHbix Bejirnvm b 
ny6jiHKaLi,Hflx h na caiiTax. 



(c) rocyflapcTBeHHbift HaynHbiii i];eHTp 
PoccHHCKofi Qe^epai^HH 

HncTHTyT (J)H3HKH BbicoKHx SHeprHftState Research Center of Russia 
Institute for High Energy Physics, 2003 



1 Motivation 



Fundamental physical constants (FPC) are the basic entities in pure and applied natural sciences 
and in technology. Thanks to efforts of many national metrology institutions, NIST^ systematists, 
and international coordination by CODATA^ we seem to have a more or less reliable procedures to 
monitor the development of the unified system of constants, their periodically adjusted numerical 
values, uncertainties, and correlations. 

We use "seem" because in spite of many national and international documents on the rules and 
standards on the statistical and experimental data presentations in the official publications^ there 
is a sharp contradiction between rules /standards and reality of scientific data exchanges in the past 
and modern scientific communication media. 

An example of such contradictions is the lack of attention from producers (NIST) and over- 
seers (CODATA) of the evaluated FPC data to the quality of the final data presentations in the 
official publications on the paper and even in electronic forms: data presented are incomplete and 
inaccurate as we will show further. 

The other example of the mentioned contradiction is the ignorance of correlations in all re- 
spectable information resources: handbooks, textbooks, monographs, reviews, and scientific soft- 
ware packages that have reprinted samples of the recommended FPC-1998 (see, for example, [5] - 

[11])- 

Since the release of FPC-1998 the scientific community obtained real access to the correlations 
of the FPC uncertainties and one can see that correlations between uncertainties of some universal 
constants are too "strong" to be ignored in high accuracy calculations. But unfortunately the 
correlation matrix presented on the NIST/CODATA site and reproduced partly in the publication 
[1] (see. Table XXV on the page 453) is non positive semi-definite in contradiction with definition 
of the correlation matrix. Looking at the NIST correlation coefficients one can see that they are 
rounded off too "tightly." Format of the numbers shows that they were rounded uniformly and 
independently (e.g. ignoring the crucial constraints that any covariance and correlation matrix 
must respect). 

^National Institute of Standards and Technology(USA) [2]. 
^Committee on Data for Science and Technology [3]. 

^We failed, however, to find any official documents standardising the procedures of rounding average values, 
standard uncertainties, and their correlations of the jointly measured or evaluated (adjusted) quantities. 
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Table 1: Sample of a few NIST/CODATA:1998 recommended constants. 



FPC name 


Simbol [units] 


Value (uncertainty) x scale 


Correlations 


Elementary charge 
Planck constant 
Electron mass 
Proton mass 


e [C] 
h [Js] 
rUe [kg] 

TTlp [kg] 


1.602176462(63) x lO"^'^ 
6.626 068 76(52) x lO'^^ 
9.109 38188(72) x 10'^^ 
1.672 62158(13) x lO"^^ 


e h rrie nip 


0.999 

0.990 0.996 

0.990 0.995 1.000 



From the textbooks on numerical calculations it is known that the rounding of the correlated 
values is a subject of special treatment (see for example [12], page 499). The average values, 
standard uncertainties, and correlation coefficients could not be rounded off independently. 

Independent rounding may lead to catastrophic changes in the connection of averages, standard 
uncertainties, and scatter ellipsoid: average values may get out of scatter ellipsoid, scatter ellipsoid 
may turn to hyperboloid after independent rounding off the correlation coefficients. 

With the lack of discussions of the rounding correlated quantities in the NIST/CODATA pub- 
lications, we interpret this as the result of independent rounding that destroy catastrophically 
the system of adjusted values. An example of such catastrophe with FPC-1998 recommended for 
"public usage" is the negative variance for the Rydberg constant calculated from the equation 
Roo = a^mcc/2h [13]. If this confusion did not caused by a misprint in signs of some correlation 
coefficient it is most probably the "inconsistency induced by rounding off" the adjusted FPC-1998. 

If we are right in our account that corruption of the true data was due to independent rounding 
then we should clarify the influence of independent rounding of the average values and dispersions. 
To test this a reproduction of the whole adjustments procedure used by NIST experts is needed. 

The main goal of this work is to reproduce the NIST results using their data and methods as 
they presented in the detailed publication [1 ] and on the NIST site and then to work out the proper 
way of correlated data presentation and exchange. 

The other goal of our work is to draw attention (once again) of the physics community to the 
problem with standardization of the statistical (experimental) data presentation in the modern 
scientific communication media. 

The rest of the paper is organized as follows. In the next section we present a few variants 
of our adjustment in comparison with corresponding NIST/CODATA results. The last section 
summarizes outputs from our exercises and our vision how to improve the situation. 

2 On the NIST FPC Adjustments Technology 

In our analyses we tried to be as close as possible to the NIST adjustments strategy. Fortunately 
more or less complete overview of the NIST: 1998 adjustment procedure including detailed pre- 
sentation of the experimental data, theoretical models and formulae, the FPC evaluation strategy 
description, and explanation of the specific aspects of the calculations were published in [1]. To test 
the traceability of the NIST results through their public information resources we have attempted 
to reproduce NIST FPC values on the basis of data, formulae and instructions from publication [J ] 
and NIST site only. 
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Table 2: Basic adjusted constants. 





Symbol 


TTPP 1 QQS Trc»liioc A/*^ 9 
rrK^-ltyuo vaiueb VO.Z 


Connections 


1 




10 973 731.568 549(83) 


_Roo = a^mcc/2h 


2 




5.485 799110(12) x 10"-* it 




3 




1.007 276 466 88(13) ?i 




4 




1.008 664 915 78(55)11 




5 


Ar{d) 


2.013 553 212 71(35) It 




u 








7 


Ar{a) 


4.001 506 1747(10) u 




8 


a 


7.297 352 533(27) x 10"^ 


a — /2eo he 


y 




— DOo.ziU Do/0(^DDj 




1 n 


Md/ Me 


— A fifi4 "^4^ ^'^71' ^0"! V 1 O"'^ 

'4:.UU'4: 0'4:tJ OO 1 yOKJ J A _LU 




]^ 


Pe- / Mp 


— fi^S 997 ^Q^d(7^ ^ 




1 9 


//' ///' 


7fi1 7kf^ 1 "^l '^T'^'^'l 
U. 1 Ul ( OU ±o±o\^oo J 




1 








14 




4.836 332 10(15) x 10^^ 




15 


h 


6.626 068 76(52) x lO"^"* Js 




16 


R 


8.314 472(15) J mol-^K-i 




1 7 


ALII \_>U. iVLt]. I 


1 nn9 n77fl'H('9s'i x in~i^m 




18 


xu(Mo Kqi) 


1.002 099 59(53) x 10"" m 


A(MoKci) = 707.831xu(Mo Kai) 


19 


A* 


1.000 015 01(90) X 10""' m 


\(WK^i) = 0.2090100A* 


20 


d22o(ILL) 


not given 




21 


d220(N) 


not given 




22 


d22o(W17) 


not given 




23 


d22o(W04) 


not given 




24 


d22o(W4.2a) 


not given 




25 


d22o(MO*4) 


not given 




26 


d22o(SHl) 


not given 




27 


d220 


1.920155 845(56) x 10"^° m 




28 


Rp 


0.907(32) X 10""' m 


From [i ], page 440 


29 


Rd 


2.153(14) X lO"""^ m 


From [ ] , page 440 



There are 107 principal input data expressed as functions of 57 adjusted constants (variables) 
via the set of observational equations. These 57 variables are subdivided in two classes: constants 
(see Table 2) and corrections (see Table 3) 

Numerical values in the Tables 2, 3 are presented in the standard concise form X{Y) x 10^ u 
(where X — average value, Y — standard uncertainty of X referred to the last digits of the quoted 
value, u — unit) and are expressed in SI. The symbols of adjusted constants are as follows: Roo 
(Rydberg constant); a (fine structure constant); h (Plank constant); mg/m^ (electron-muon mass 
ratio); Ar{X) ( mass of the particle X in atomic units), where X denotes symbol of the particle, 
such as electron or alpha-particle; fJ.x/fJ'Y {X, Y magnetic moment ratio); R — molar gas constant; 
xu(XKai) — x-ray unit of the X atom; d22o(X) — {220} lattice spacing of the different (X) silicon 
mono-crystal (d22o); Rp and Rd are the bound-state proton and deuteron rms charge radii. 

There are also 28 adjusted variables (see Table 3) that are not fundamental at all, but were 
introduced in order to decrease the theoretical uncertainty of the several observational equations. 
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Table 3: Basic adjusted corrections. 





INibi bymbol 


INlb i /L-(JUAlA:19yo Value 


Comments 


30 


Sh{1,0,1) 


not given 




31 


5^(2,0,1) 


not given 




32 


(5//(3,0, 1) 


not given 




33 


5j/(4,0,l) 


not given 




34 


5_f/(6,0, 1) 


not given 




35 


5^(8,0,1) 


not given 




36 




not given 




37 




not given 




38 


fe(2,l,3) 


not given 




39 


fo(4, 1,3) 


not given 




40 


i5_f/(8,2,3) 


not given 




41 


5//(12,2,3) 


not given 




42 


5//(4,2,5) 


not given 




43 


i5_f/(6,2,5) 


not given 




44 


(5_f/(8,2,5) 


not given 




45 


5//(12,2,5) 


not given 




46 


(5£,(1,0,1) 


not given 




47 


(5i,(2,0,l) 


not given 




A Q 


1)0(4, 0, Ij 


not given 




49 


(5d(8,0,1) 


not given 




50 


(5d(8,2,3) 


not given 




51 


fo(12,2,3) 


not given 




52 


5d(4,2,5) 


not given 




53 


(5d(8,2,5) 


not given 




54 


fo(12,2,5) 


not given 




55 


5e 


0.1(1.1) X 10-12 


From [1], page 457 


56 




not given 




57 




0.0(6.4) X 10-1" 


From [1], page 457 



They are: 

6TSs{n, L, 2 • j) — additive correction to riLj energy level of the hydrogen (N=H) and deuterium 
(N=D); 

Se, Sfi — additive correction to electron and muon magnetic moment anomaly, 
^Mu — additive correction to hyperfine splitting of the muonium basic bound state. 

Unfortunately the set of final values released so far by NIST/CODATA:1998 does not contain 
estimates for large part of the 57 adjusted variables (marked as "not given" in the Tables 2, 3 
in cases when we failed to find corresponding output value on the site CODATA Fundamental 
Physical Constants. Version 3.2 Release date: 1 October 2003. [1]). Hence our comparison will be 
incomplete to this extent. 
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The other problem that makes the straightforward reproduction of the NIST estimates impos- 
sible is the corrupted presentation of the input correlation sub-matrix (see [i], page 434, Table 
XIV. A. 2) of the data related to the Rydberg constant. It is non positive definite. We interpret this 
confusion as the result of unjustified independent rounding of correlation coefficients (motivated 
only to be convenient for publication on the paper). Fortunately the correlation sub-matrix for the 
other data presented in [1] (see page 436, Table XIV. B. 2 ) is positive definite. 

In spite of the incomplete presentation of the adjustment results and corrupted input correlation 
data we decided to clarify to what extent the ignorance of the input correlations in part or totally 
will modify the output constants, supposing that average values of the constants are correct as well 
as their standard uncertainties. 

The strategy of comparisons is as follows. First of all we have convinced that the NIST method 
used to obtain values of adjustable variables is indeed the method to find stationary points of the 
"linearized" x^- If the input covariance matrix is positive definite then the obtained solution will be 
a minimum and adjustment will be stable if we turn lucky to get into vicinity of a global minimum. 
If the input covariance matrix is non positive definite then the task of finding stationary points 
could be solved for non-degenerate weight matrices but in this case the task has no connection with 
least squares method of constants estimation. 

Nevertheless it is interesting to compare values of adjusted variables that can be obtained by 
NIST method in our PAREVAL package (from the NIST starting point in the adjusted variables 
space and with rounded by NIST input covariance matrix [1]) with supposedly true NIST:1998 
values. 



2.1 Calculations with input correlator published by NIST 

Corresponding results are presented in Table 4, where in the fifth column the normalized difference 
of values of third (X^) and fourth (X4) columns is defined as follows 



In our adjustment we use close to the NIST starting point (used for finding zero of the gradient 
of in step-by-step method) in the adjusted variables space and "weights" were constructed from 
rounded by NIST input "correlation" matrix [1]. 

As it can be seen from the fifth column of the Table 4, the shifts in average values are in 
general well inside the ranges defined by the quadratically combined "uncertainties". In all tables, 
as a rule, we save in average values and corresponding uncertainties one more digit to show that 
some values are not reproduced exactly when rounding independently even if they are well inside 
the uncertainties (|A| < 0.1). 
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Table 4: Comparison of the true NIST/CODATA:1998 recommended values (third column) with corre- 
sponding values (fourth column) that have been obtained at IHEP by the NIST method in our PAREVAL 
package. 





FPC 


NIST: 1998 Value 


IHEP: 2003 Value 






Symbol 


NIST true correlator 


NIST published correlator 


A 


1 




1.097373156 8549(83) x 10'' 


1.097373156 854 7(83) x lO'^ 


0.0153 


2 


Ar.(e) 


5.485 799110(12) x 10"^ 


5.485 7991097(116) x IQ-"^ 


0.0171 


3 


Ar.(v\ 


1.007276 466 88(13) 


1.007 276 466 883(132) 


-0.0153 


4 


Ar{n) 


1 008 664915 78(55) 


1 008 664 915 784(547) 


00501 


5 


Ar{d) 


2.013 553 212 71(35) 


2.013 553 212 706(344) 


0.00833 


6 


Ar{h) 


3 014932 23469(86) 


3 014932 234691(860) 


-0.001 


7 


Ar{a) 


4.0015061747(10) 


4.00150617469(100) 


0.00456 


8 


Ct 


7.297352 533(27) x 10~^ 


7.297352 5335(265) x IQ-^ 


-0.0132 


9 


n„— 1 n„ 
H'e 1 H'P 


-6.582106 875(66) x 10^ 


-6.582106 8753(659) x 10^ 


0.00375 


10 


UrI / U~ 
r'CLI r^e 


-4.664345 537(50) x 10"^ 


-4.664 345 5371(500) x IQ-"^ 


0.0009 


11 


Li — 1 u! 


-6.582 275 954(71) x 10^ 


-6.582 275 9549(717) x 10^ 


0.00938 


12 


Hhl Hp 


-7.617861313(33) x 10"^ 


-7.6178613130(330) x IQ-i 


6.64 X 10-^ 


13 


n,^ 1 ti! 
Hn/ Hp 


-6 849 9694(16) x 10^^ 


-6 849 96940(160) x 10"! 


6.16 X IQ-i^ 


14 




4 836 332 1 on 5^1 x 10"-^ 


4 836 33? 107(144') x 10-^ 


—0 0360 


15 


h 


6 626 068 76(52) x 10""^^ 


6 626 068 756(522") x lO"'^"' 


00530 


16 


R 


8.314472(15) 


8.3144724(147) 


-0.0214 


17 


xu(Cu Kai) 


1.002 07703(28) x lO"^^ 


1.002 077021(287) x lO-^^ 


0.0212 


1 Q 


xul^lVlo J\ai ) 


i.uuz uyy oy(t)o) x lu 


1 nno nnn c^oQ 1 ^/ in~13 

i.uuz (jyy oyo(oioj x lu 


— U.UU4Di 


19 


A* 


1.000 015 01(90) X 10-10 


1.000 015 010(901) X 10-1° 


-0.000017 


20 


d22o(ILL) 


not given 


1.920155 8160(558) x 10-i° 




21 


d220(N) 


not given 


1.920155 8191(508) x 10-i° 




22 


d22o(W17) 


not given 


1.920155 8411(484) x 10-i° 




23 


d220(W04) 


not given 


1.920155 8103(520) x 10-1° 




24 


d22o(W4.2a) 


not given 


1.920155 7995(496) x 10-i° 




25 


d220(MO*4) 


not given 


1.920155 6075(439) x 10-i° 




26 


d22o(SHl) 


not given 


1.920155 7624(463) x 10-i° 




27 


d220 


1.920155 845(56) x 10"^° 


1.920155 8391(561) x 10-i° 


0.0749 


28 


Rp 


0.907(32) X 10-15 


0.9066(329) x lO-i^ 




29 


Rd 


2.153(14) X 10-15 


2.1528(137) X 10-15 








x'^/ndf = 0.90 [1], page 446 


^Indj = 0.90 
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As we have no NIST/CODATA:1998 final values of the adjusted corrections we compare our results 
with the input values which we use (follow NIST) as observational equation like S = 0.0 . . . 0{u{5)) 
(see Table 5). We see that shifts in average values for all corrections are also well inside the 
intervals defined by the input estimates of theoretical systematic uncertainties. 



Table 5: Comparison of the input NIST/CODATA:1998 values (third column) with corresponding values 
(fourth column) that obtained at IHEP by the NIST method in our PAREVAL package. 





FPC 
Symbol 


NIST:1998 input value 


IHEP: 2003 value 
NIST published correlator 


30 


(5h(1,0, 1) 


0.0000(9.0000) X 10" 


0.009(8.974) X 10" 


31 


fe(2,0,l) 


0.0000(1.1000) X 10'' 


0.007(1.097) X 10" 


32 


fe(3,0,l) 


0.000(3.300) X 10^ 


0.003(3.291) X 10^ 


33 


(5h(4,0,1) 


0.000(1.400) X 10^ 


0.009(1.396) X 10^ 


34 


ShIg, 0,1) 


0.00(4.20) X 10^ 


0.00(4.19) X 10^ 


35 


(^H 8,0,1) 


0.00(1.80) X 10^ 


0.00(1.79) X 10^ 


36 


'5h(2,1,1) 


0.000(1.100) X 10^ 


-0.04(1.09) X 10^ 


37 


-^^(4,1,1) 


0.00(1.40) X 10^ 


-0.05(1.39) X 10^ 


38 


(5^(2,1,3) 


0.000(1.100) X 10^ 


0.01(1.10) X 10^ 


3Q 




n onn 40"! x i o^ 


n C\9h 391 X 1 0^ 


40 


(5h(8,2,3) 


0.0(1.7) X IQi 


-0.00(1.70) X 10^ 


41 


5h(12,2,3) 


0.0(5.0) 


-0.007(5.000) 


42 


fe(4,2,5) 


0.00(1.40) X 10^ 


0.01(1.40) X 10^ 


43 


fe(6,2,5) 


0.0(4.0) X 10^ 


0.03(4.00) X 10^ 


44 


<5h(8,2,5) 


0.0(1.7) X IQi 


0.10(1.70) X IQi 


45 


<5h(12,2,5) 


0.0(5.0) 


0.04(5.00) 


46 




0.0000(8.9000) X 10" 


0.005(8.877) X 10" 


47 


-5^(2,0,1) 


0.0000(1.1000) X 10" 


0.001(1.097) X 10" 


48 


5c (4, 0,1) 


0.000(1.400) X 10^ 


0.001(1.396) X 10^ 


49 


SniS, 0,1) 


0.00(1.70) X 10^ 


0.00(1.70) X 10^ 


50 


(5z5(8,2,3) 


0.0(1.1) X IQi 


-0.00(1.10) X 10^ 


51 


(5/3(12,2,3) 


0.0(3.4) 


-0.00(3.40) 


52 


fo(4,2,5) 


0.0(9.2) X 10^ 


0.06(9.20) X 10^ 


53 


Sd{8,2,5) 


0.0(1.1) X IQi 


0.008(1.100) X lOi 


54 


fo(12,2,5) 


0.0(3.4) 


0.02(3.40) 


55 


6e 


0.00(1.1) X 10-12 


0.01(1.10) X 10-12 


56 




0.00(1.20) X 10^ 


0.01(1.17) X 102 


57 




0.0(6.4) X 10-1" 


0.00(6.40) X 10-1° 








X^/ndf = 0.90 



7 



2.2 Calculations with identity input correlation matrix 

We learn from Tables 4 and 5 that rounded input correlation matrix leads to slightly shifted values of 
the basic constants and corrections that are well inside the uncertainties. To clarify the importance 
of correlations further and to test our package we produce adjustment ignoring input correlations. 



Table 6: Comparison of the true NIST/CODATA:1998 recommended values (third column) with corre- 
sponding values (fourth column) that obtained at IHEP by the NIST method in our PAREVAL package. In 
this adjustment we use the same starting point as in previous run and weights were constructed completely 
ignoring input correlations. We see (from the fifth column) that shifts in average values are in general inside 
the quadratically combined "uncertainties" . 





FPC 


NIST: 1998 Value 


IHEP: 2003 Value 






Symbol 


NIST true correlator 


Identity matrix correlator 


A 


1 


R 


1.097 373156 8549(83) x lO'^ 


1.097 373156 8545(103) x lO'^ 


0.0314 


2 


Ar{e) 


5.485 799110(12) x 10^* 


5.485 7991099(116) x 10"" 


0.00317 


3 




1.007 276 466 88(13) 


1.007 276 466 883(132) 


-0.0138 


4 


Ar{n) 


1.008 664 915 78(55) 


1.008 664 915 774(556) 


0.00731 


5 


Ar{d) 


2.013 553 212 71(35) 


2.013 553 212 688(360) 


0.0439 


r- 
D 


Ar{ri) 


3.014 \)6l z34 D9(86j 


3.014 y3z 234 d9U(86U ) 


U.0004U2 


7 


Ar{a) 


4.0015061747(10) 


4.001 506 17469(100) 


0.00459 


8 


a 


7.297 352 533(27) x 10"^ 


7.297 352 5349(266) x lO"'' 


-0.0500 


9 




-6.582 106 875(66) x 10^ 


-6.582 106 8754(659) x 10^ 


0.00432 


10 




-4.664 345 537(50) x lO"* 


-4.664 345 5371(500) x 10"'' 


0.000901 


11 




-6.582 275 954(71) x 10^ 


-6.582 275 9550(717) x 10^ 


0.00973 


12 




-7.617 861313(33) x 10'^ 


-7.617 8613130(330) x 10"^ 


6.64 X lO"'' 


13 


-6.849 9694(16) x 10"^ 


-6.849 96940(160) x 10"^ 


6.16 X 10"^^ 


14 


me/m^ 


4.836 332 10(15) x 10"^ 


4.836 332 103(142) x 10"=* 


-0.0124 


15 


h 


6.626 068 76(52) x 10"^-* 


6.626 068 756(522) x 10"=''' 


0.00523 


16 


R 


8.314 472(15) 


8.314 4724(147) 


-0.0214 


17 


xu(Cu Koi) 


1.002 07703(28) x 10"" 


1.002 077018(288) x 10"^=* 


0.0306 


18 


xu(MoKai) 


1.002 099 59(53) x 10"^^ 


1.002 099 596(516) x 10"^^ 


-0.00803 


19 


A* 


1.000 015 01(90) X 10"^" 


1.000 015 013(901) X 10""' 


-0.00200 


20 


d22o(ILL) 


not given 


1.920 155 8093(517) x 10"^° 




21 


d220(N) 


not given 


1.920 155 8239(635) x 10"^° 




22 


d22o(W17) 


not given 


1.920 155 8380(545) x 10"^° 




23 


d22o(W04) 


not given 


1.920 155 8102(460) x 10"^" 




24 


d22o(W4.2a) 


not given 


1.920 155 7910(563) x 10"^° 




25 


d22o(MO*4) 


not given 


1.920 155 5957(457) x 10"^° 




26 


d22o(SHl) 


not given 


1.920 155 7631(509) x 10"^° 




27 


d220 


1.920 155 845(56) x 10"" 


1.920 155 8390(506) x 10"^° 


0.0800 


28 


Rp 


0.907(32) x 10"^'^ 


0.9087(261) X 10"^^ 




29 


Rd 


2.153(14) X 10"^^ 


2.1519(115) X 10"^^ 








X^/ndf = 0.m [I], page 446 


X^/ndf = 0.84 
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Table 7: Comparison of the obtained values of additive corrections with input estimates of the 
theoretical systematic uncertainties. 





r r U 
bymboi 


iNioi:iyyo input value 


iJlr!jr^:zUUo value 
Identity matrix correlator 


30 


(5^(1,0,1) 


0.0000(9.0000) X 10-^ 


-0.06(6.40) X 10^ 


31 


6h{2,0,1) 


0.0000(1.1000) X 10^ 


-0.07(7.96) X 10^ 


32 


Sh{3,0,1) 


0.000(3.300) X 10^ 


0.01(3.29) X 10^ 


33 


(5^(4,0,1) 


0.000(1.400) X 10^ 


0.01(1.39) X 10^ 


34 


5i?(6,0,l) 


0.00(4.20) X 10^ 


0.02(4.20) X 10^ 


35 


Sh{8,0,1) 


0.00(1.80) X 10^ 


-0.01(1.80) X 102 


36 


Sh{2,1,1) 


0.000(1.100) X 10^ 


-0.05(1.10) X 10^ 


37 




0.00(1.40) X 10^ 


0.01(1.40) X 10^ 


38 


fe(2,l,3) 


0.000(1.100) X 10^ 


0.03(1.10) X 10^ 


39 


<5i?(4,l,3) 


0.00(1.40) X 10^ 


-0.02(1.40) X 10^ 


40 




0.0(1.7) X 10^ 


0.00(1.70) X 10^ 


41 


fe(12,2,3) 


0.0(5.0) 


-0.00(5.00) 


42 


(5/^(4,2,5) 


0.00(1.40) X 10^ 


0.005(1.400) X 102 


43 




0.0(4.0) X 10^ 


-0.009(4.000) X IQi 


44 


<^//(8,2,5) 


0.0(1.7) X 10^ 


0.000(1.700) X IQi 


45 


Sh{12,2,5) 


0.0(5.0) 


-0.00(5.00) 


46 


<5i5(l,0,l) 


0.0000(8.9000) X 10"^ 


0.05(6.50) X 10"^ 


47 


Sd{2,0,1) 


0.0000(1.1000) X 10^ 


-0.04(8.07) X 10^ 


48 




0.000(1.400) X 10^ 


-0.08(1.40) X 10^ 


49 


6d{8,0,1) 


0.00(1.70) X 10^ 


0.02(1.70) X 10^ 


50 


fe(8,2,3) 


0.0(1.1) X 10^ 


-0.00(1.10) X 10^ 


51 


<^d(12,2,3) 


0.0(3.4) 


-0.00(3.40) 


52 


<^d(4,2,5) 


0.0(9.2) X 10^ 


0.00(9.20) X 10^ 


53 


6d{8,2,5) 


0.0(1.1) X 10^ 


0.00(1.10) X 10^ 


54 


Sd{12,2,5) 


0.0(3.4) 


-0.00(3.40) 


55 


Se 


0.0(1.1) X 10-12 


0.08(1.10) X 10-12 


56 




0.00(1.20) X 10^ 


0.01(1.17) X 10^ 


57 




0.0(6.4) X 10-1° 


0.00(6.40) X 10-1° 








x''/ndf = 0.84 
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Now we can compare values of constants and corrections obtained with rounded correlator and 
without correlations (see Tables 8, 9). 



Table 8: Comparison of constant estimates obtained with NIST published input correlations and 
with identity matrix as input correlation matrix 





FPC 


IHEP:2003 value 


IHEP:2003 Value 






Symbol 


NIST published correlator 


Identity matrix correlator 


A 


1 




1.097 373 156 855(8) x 10 


1.097 373 156 854(10) x 10 


0.018 


2 


Up 


9.066(329) X 10 


9.087(261) x 10 


—0.048 


6 


Rd 


2.1528(137) X 10 


2.1519(115) X 10 


0.048 


A 

4 


Ar[e) 


5.485 799 1097(116) x 10 


5.485 799 1099(116) X 10 


n m A 

—0.014 


5 


A / \ 

Ar[p) 


1.007 276 466 883(132) 


1 r\r\''7 o*^/^ A r* n 000/1 ooN 

1.007 276 466 883(132) 


0.0015 


6 


A / \ 

Ar[n) 


1.008 664 915 784(547) 


1 r\r\c> CCA (\'\ v A I v c r'\ 

1.008 664 915 774(566) 


0.012 


''7 

7 


Ar{a) 


2.013 553 212 706(344) 


c\ n\-\ r'r'o 010 /^oo/o/^/~\\ 

2.013 553 212 688(360) 


0.036 


o 





3.014 932 234 691(860) 


3.014 932 234 690(860) 


0.0014 


9 


A / \ 


4.001 506 174 69(100) 


4.001 506 174 69(100) 


0.000 33 


10 


a 


7.297 352 5335(265) x 10 


7.297 352 5349(266) x 10 


—0.037 


11 




—6.582 106 8753(659) x 10 


—6.582 106 8754(659) x 10 


0.0058 


12 


/ — 


—4.664 345 5371(500) x 10 


A CCA O/ir' r'o'7i/r'r\/~\N ir\ — 4 

—4.664 345 5371(500) x 10 





13 


1 ,,t 

IJ'e-/fJ'p 


—6.582 275 9549(717) x 10 


—6.582 275 9550(717) x 10 


0.0035 


14 


/ / f 


—7.617 861 3130(330) x 10 


^ c '{ ^ 0/^1 010^\/oo/~^^ ir\- — 1 

—".Gil 861 3130(330) x 10 





iO 


II III' 

^l-nlU'p 


fi Q/1Q QfiQ /in^'i w 1 n^i 
— D.o4y yoy 4U(^iDUj x lu 


— D.o4y yoy 4U(^iDUj x lu 


n 
u 


16 




4.836 332 107(144) x 10"^ 


4.836 332 103(142) x 10^^ 


0.024 


17 


h 


6.626 068 756(522) x 10^34 


6.626 068 756(522) x 10^34 


-0.000 64 


18 


R 


8.314 4724(147) 


8.314 4724(147) 





19 


xu(Cu Kai) 


1.002 077 021(287) x lO'^^ 


1.002 077 018(288) x lO^^'^ 


0.0093 


20 


xu(Mo Ktti) 


1.002 099 593(516) x lO'^^ 


1.002 099 596(516) x lO^^^ 


-0.0035 


21 


A* 


1.000 015 010(901) X 10^1° 


1.000 015 013(901) X 10^1° 


-0.0020 


22 


d22o(ILL) 


1.920 155 8160(558) x 10"!° 


1.920 155 8093(517) x 10"!° 


0.087 


23 


d220(N) 


1.920 155 8191(508) x 10"^° 


1.920 155 8239(635) x 10"!° 


-0.060 


24 


d22o(W17) 


1.920 155 8411(484) x 10"^° 


1.920 155 8380(545) x 10"!° 


0.043 


25 


d220(W04) 


1.920 155 8103(520) x 10"^° 


1.920 155 8102(460) x 10^^° 


0.0015 


26 


d22o(W4.2a) 


1.920 155 7995(496) x 10"^° 


1.920 155 7910(563) x 10^^° 


0.011 


27 


d220(MO*4) 


1.920 155 6075(439) x 10"^° 


1.920 155 5957(457) x 10"!° 


0.019 


28 


d22o(SHl) 


1.920 155 7624(463) x 10"^° 


1.920 155 7631(509) x 10"!° 


-0.0096 


29 


d220 


1.920 155 8391(561) x 10"^° 


1.920 155 8390(506) x 10"^° 


0.0013 






^Indj = 0.90 


X^/ndf = 0.84 
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Table 9: Comparison of correction estimates obtained with NIST published input correlations and 
with identity matrix as input correlation matrix. 





FPC 
Symbol 


IHEP:2003 value 
NIST published correlator 


IHEP:2003 Value 
Identity matrix correlator 


A 


30 


6h(1,0, 1) 


0.009(8.974) X 10^ 


-0.062(6.40) X 10^ 


0.057 


31 


5ui2 1) 


007(1 0974) X 10^ 


-0 067(7 96) X 10^ 

\J * \J \J 1 \ 1 ■ t/ \J } /^N -L KJ 


0.050 


32 


5i/(3,0, 1) 


0.003(3.291) X 10^ 


0.012(3.29) X 10^ 


-0.024 


33 


5^(4,0, 1) 


0.009(1.3963) X 10^ 


0.014(1.39) X 10^ 


-0.069 


34 


(5w('6 11 


0f4 19) X 10^ 


0.02(4.20) X 10^ 


-0.0027 


35 


5h(8,0, 1) 


0(1.79) X 10^ 


-0.01(1.80) X 10^ 


0.0056 


36 


5h{2, 1, 1) 


-0.04(1.09) X 10^ 


-0.05(1.10) X 10^ 


0.0069 


37 


5ff(4, 1, 1) 


-0.05(1.39) X 10^ 


0.01(1.40) X 10^ 


-0.032 


38 


(Jw(2, 1,3) 


0.01(1.10) X 10^ 


0.03(1.10) X 10^ 


-0.011 


39 


5f7('4, 1,3) 


0.02(1.39) X 10^ 


-0.02(1.40) X 10^ 


0.019 


40 


5^(8,2,3) 


-0(1.70) X 10^ 


0(1.70) X 10^ 


-0.0013 


41 


5w('12,2,3) 


-0.007(5.000) 


-0(5.00) 


-0.0060 


42 


(5w('4, 2,5) 


0.01(1.400) X 10^ 


0.005(1.400) X 10^ 


0.0024 


43 


5wf6 2 5) 


03f4 00) X 10^ 

KJ • KJ y \ i a yj KJ J / V _L. KJ 


-0 009f4 000) X 10^ 

KJ • KJ \J \J I -1- • KJ \J KJ } X KJ 


0065 

KJ • KJ \J KJ \J 


44 


5jy(8,2,5) 


0.01(1.70) X 10^ 


0(1.70) X 10^ 


0.0030 


45 


5j/(12,2,5) 


0.04(5.00) 


-0(5.00) 


0.0053 






n nn/iQ/'s R7f\7\ v in^ 
u.uu4y i^o.o ( ( j X iu 


U.U4D\^D.0Uj X iU 


n D/i 1 


47 


5d(2,0, 1) 


0.0011(1.0970) X 10^ 


-0.04(8.07) X 10^ 


0.0034 


48 


fo(4,0,l) 


0.001(1.396) X 10^ 


-0.08(1.40) X 10^ 


0.043 


49 


fo(8,0, 1) 


0(1.70) X 10^ 


0.02(1.70) X 10^ 


-0.0076 


50 


fe(8,2,3) 


-0(1.10) X 10^ 


0(1.10) X 10^ 


-0.0012 


51 


fe(12,2,3) 


-0(3.40) 


-0(3.40) 


-0.0080 


52 


fo(4,2,5) 


0.06(9.20) X 10^ 


0(9.20) X 10^ 


0.0045 


53 


fe(8,2,5) 


0.008(1.100) X 10^ 


0(1.10) X 10^ 


0.0039 


54 


fe(12,2,5) 


0.02(3.40) 


-0(3.40) 


0.0052 


55 




0.010(1.10) X 10~i^ 


0.08(1.10) X 10~i2 


0.0091 


56 




0.013(1.17) X 10^ 


0.012(1.17) X 10^ 


0.0058 


57 




0(6.40) X 10^1° 


0(6.40) X 10-1° 


0.000 20 






x''/ndf = 0.90 


^Indj = 0.84 
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Finally we have checked if the procedure to find minimum of the corresponding without 
linearization of the observational equations used gives the same result. We have used the built-in 
Mathematica module to find minimum starting from the point in the adjusted variables space where 
the values of all 29 constants were taken as recommended by NIST/CODATA:1998 and the values 
of the rest 28 (5-corrections were taken as 0.001 x u{6theor)- 

Corresponding value of /ndf = 0.84, i.e. the same as obtained by using NIST procedure with 
linearization. Average values of the constants and corrections obtained by two different methods and 
different programs are practically the same (the maximal normalized difference is |A|maa; ~ 10~^^). 

Our results lead us to the conclusion that NIST experts have used input correlation matrix 
close to that of presented in their published report. The fact that output correlation matrix is 
also non positive semi-definite does not allow one to exclude possibility that this non positive 
semi-definiteness could be induced by the non positive definiteness of the input matrix. This point 
remains to be clarified. 

Let us compare a sample of our output correlation coefficients with corresponding NIST values. 



Table 10: The values in bold are extracted from the NIST site whereas the values placed under 
main diagonal are our values obtained with no input correlations. 





Roo 


a 


me/nifj, 


h 




1.00 


-0.020 


0.004 


-0.000 


a 


-0.0112 


1.00 


-0.233 


0.002 




0.00235 


-0.236 


1.00 


-0.000 


h 


-0.0000195 


0.00174 


-0.000410 


1.00 



As shown in the Table 10 some of our correlation coefficients differ significantly from those of 
pubhshed by NIST. 

Now we proceed to the general comments on the practice in the scientific literature and on the 
sites, where the correlated estimates of the random quantities jointly measured or evaluated are 
presented. 

Often authors of the original papers appear to ignore correlations at all or to present them in an 
incomplete manner. So, it is hard to understand what type of uncertainties the quoted correlation 
matrix is referred to: statistical, systematic or total. It is dangerous (or even inadmissable) to use 
such incomplete data in further analyses and especially in the theory tests. Sometimes this incom- 
pleteness is caused by the too firm editors and publishers requirements. On the other hand there 
are also experimental mistakes (see examples of such situation in nuclear physics and technology 
[17]). 

2.3 On the rounding off the correlated estimates 

As it is well known the covariance matrix for the jointly estimated statistical quantities is by 
definition a positive semi-definite Hermitian matrix. It has real and nonnegative eigenvalues. The 
non degenerate correlation matrix is positive definite by definition. 
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Unfortunately some authors publish the "correlation matrix" with no final check up this crucial 
property of the correlation matrix. In majority cases it happens under the pressure of the limited 
publication space. So, authors are forced to present rounded correlation coefficients making this in 
an inadmissable manner. The rounding off the correlation coefficients are produced independently, 
saving symmetry of the matrix but ignoring such crucial properties as positive definiteness and 
positive semi-definiteness. 

To be specific we quote a few examples from different subject fields. The first most striking 
example is yet discussed concerning the NIST publications on the adjusted fundamental constants, 
including the NIST site. The published version of the input correlation sub-matrix used to construct 
the weight matrix in their version of least squares method (LSM) is non positive definite, (see 
[1], page 434, Table XIV. A. 2), it has two negative eigenvalues. Also the published version of the 
correlation sub-matrix between uncertainties in the recommended values of a sample of fundamental 
constants (see [1], page 453, Table XXV) is non positive semi-definite. We have convinced that it is 
because of poor accuracy of the presentation caused by unjustified uniform independent rounding 
of the correlation coefficients. Unfortunately the same way of presentation is used by NIST and 
approved by CODATA on their sites. 

The other example is the publication of the CLEO collaboration on the high precision mea- 
surements of the r-lepton decay branching ratios [II]. The final version of the correlation matrix 
presented in the Erratum is as follows: 

Table 11: Correlation coefficients between branching fraction measurements of some r-lepton decays 
from [14] . 



r 


Be 




Bh 


B^,/Be 


Bh/Be 


Be 


1.00 


0.50 


0.48 


-0.42 


-0.39 






1.00 


0.50 


0.58 


0.08 


Bh 






1.00 


0.07 


0.63 


Bfj./Be 








1.00 


0.45 


Bh/Be 










1.00 



The corresponding eigenvalues are as follows: 

2.173 46, 1.78187, 1.054 97, -0.00749153, -0.002 8034 

This confusion could be due to improper rounding, but we failed to show this by playing with num- 
bers (un-rounding) . The problem seems to be deeper and hence the CLEO data are questionable. 
These should be used with great caution in theory tests and in derivations of "world averaged" 
r-lepton branching ratios. 

As we already mentioned, the proper rounding procedure for the jointly measured or estimated 
quantities (average vector components, corresponding vector of their standard uncertainties, and 
correlation matrix) is the subject of special treatment. So it will be presented elsewhere if we will 
not found relevant papers published. 

In the next section we construct a simple but important estimate of the threshold accuracy 
of the correlation coefficients that should not be violated while uniform independent rounding of 
correlation matrix elements. 
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3 On the numerical presentations of correlated quantities in com- 
puter readable files and in publications 

Here we derive a simple sufficient estimate on the accuracy of a safely independent and uniform 
rounding the correlation matrix elements off. 

Let Aij be the n x n correlation matrix. It is real, symmetric, positive definite, and has matrix 
elements bounded as follows 

An = 1 for all i = I, . . . ,n and \Ai^j\ < 1.0. 

Let Bij be the "rounder" matrix, such that if it is added to the matrix Aij the obtained matrix 
Gij = Aij + Bij will be real, symmetric, positive definite and all < 1 are decimal numbers 

with k digits after the decimal point. 

It is easy to see that matrix Bij should have the following properties: 



Bii = for all i = 1, . . . ,n and \Bi^j\ < 5.0 x 10 



Let further ai < ■ ■ ■ < Un, Pi < ■ ■ ■ < Pn, and 71 < • • • < 7n be the ordered sets of eigenvalues of 
the matrices Aij , Bij , and Gij correspondingly. Then from the Weil's theorem for any I = 1, . . . ,n 
we have the following inequalities [15], [16]: 

a« + /?! < 7/ < a« + Pn- 

From the Gershgorin's theorem on the distributions of the eigenvalues of the Hermitian matrices 
[15] it follows that 

> _ 1) . 5 . lo-e^+i) = -^^^ ■ 10^^ 

and hence to have the matrix Gij as positive semi definite matrix it is sufficient to demand 

(.n-1) u 
< ai - ^ ^ ' ■ 1^-^ < 71. 

From the left inequality we have the final estimate for the threshold accuracy index for safely 
uniform independent rounding (SUIR) of the positive definite correlation matrix Aij with minimal 
eigenvalue Umin 

. n-1 

logio 



1^ ^ ^SUIR 



NOTE. According to the Weil's theorem any uniform rounding the off-diagonal matrix elements 
of the positive semi-definite covariance matrix is forbidden. 

Indeed, as rounder matrix is traceless Hermitian matrix, it obliged to have the negative minimal 
eigenvalue. Furthermore from the left inequality of the Weil's theorem statement it follows that 
any rounding could lead to the matrix with negative minimal eigenvalue. 

This note shows that the special rounding strategy should be developed^ for such covariance 
matrices as well as for the badly conditioned covariance matrices. 



*We realise that, most probably, such a strategy was developed already somewhere, but unfortunately is deeply 
hidden in the national and international metrology instructions. Some relevant information see, for example, in the 
review [17], where the analogous concerns are expressed. 
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Now we can make a further comments on the values of the FPC correlation coefficients. In 
the Table 10, to simplify comparison with NIST data, our output correlation coefficients for basic 
FPC are presented with independent rounding to three significant digits. In fact we have minimal 
eigenvalue for our output 57 x 57 correlation matrix 

amin = 1.290 020 86151811 x 10"^ 

and from the expression (1) for the critical accuracy it follows that at least 8 digits after decimal 
point should he saved, when rounding uniformly and independently, to preserve positive definiteness 
of the output correlation matrix for basic sample of FPC. 



Table 12: Comparison of a few CODATA:1998 and CODATA:2002 recommended constants. 



CODATA:1998 


Simbol [units] 


Value (uncertainty) x scale 


Correlations 


Elementary charge 
Planck constant 
Electron mass 
Proton mass 


e [C] 
h [Js] 

nie [kg] 
rrip [kg] 


1.602176 462(63) x 10"^^ 
6.626 068 76(52) x lO'^^ 
9.109 38188(72) x 10~3i 
1.672 62158(13) x lO^^^ 


e h rUe rUp 


0.999 

0.990 0.996 

0.990 0.995 1.000 


CODATA:2002 


Simbol [units] 


Value (uncertainty) x scale 


Correlations 


Elementary charge 

Planck constant 
Electron mass 
Proton mass 


e [C] 

h [Js] 
nie [kg] 
rrip [kg] 


1.602176 53(14) x 10^^^ 

6.626 0693(11) x 10"^'' 
9.109 3826(16) x lO"^^ 
1.672 62171(29) x 10^27 


e h rrie rUp 


1.000 

0.998 0.999 

0.998 0.999 1.000 



Eigenvalues of these correlation sub-matrices are as follows: 

COD AT A : 1998 {3.985, 0.0150769, 0.00536526, -0.000617335}; 

CODATA : 2002 {3.997, 0.00315831, -0.000158432, -2.83681 x IQ-^^}. 

Both matrix are non positive semi-definite, the 2002 sub-matrix is degenerate as it is seen from the 
table above. 

Another concern is the accuracy of data presentation on average values and standard uncertain- 
ties. As a rule 2002-uncertainties are more than two times larger than corresponding 1998-numbers 
(see the table above). 

It is to some extent unexpected as the NIST bibliography database on FPC contains 528 addi- 
tional to 1998 database entries classified as "experimental" and "original research" dated between 
1999 and 2002 inclusively. Unfortunately the reasons of these enlargements of the standard uncer- 
tainties compared to the V3.2(1998)-release did not commented in the V4.0(2002)~release notifi- 
cation. 



15 



4 Conclusions 



In this section we summarize the main results obtained and the discussions presented above. 

• The Mathematica package PAREVAL was created and apphed for adjustments of FPC. It 
reproduces the NIST/CODATA adjustment technology from their data and methods (as of 
1998). It is composed of a few modules to maintain: the library of theoretical models, the 
experimental data compilation; to perform FPC evaluations, and for results presentation. 
The detailed package description as well as the address for access will be presented elsewhere 
[18]. 

• With the help of PAREVAL it is shown, that the "CODATA recommended values of the 
fundamental physical constants: 1998" V3.2 and "CODATA recommended values of the fun- 
damental physical constants: 2002" V4.0 are questionable in, at least, the values of published 
correlation coefficients released in the NIST/CODATA sites. Most probably data were cor- 
rupted by unjustified rounding up the output values. 

It is argued that the released so far correlation coefficients are useful only to show the sizes 
of correlations but should not be used in the real calculations of high precision observables. 
It will be extremely useful if the released for the first time ASCII file be accompanied with 
easy computer readable files with the compact standardized names, units, average values, 
standard uncertainties and correlations presented as accurate as possible, without unjustified 
rounding up. 

• The simple estimate for the threshold accuracy sufficient for safely uniform independent 
rounding of the positive definite correlation matrix is constructed. This estimate can be used 
to trigger the data corrupted by unjustified rounding. 

• It is argued that it is an urgent need to create a common standard strategy of rounding 
interrelated quantities and common standard data structures to store and exchange the cor- 
related data in the computation media. These standards should be freely available for science, 
education and technology practitioners. 

It seems that it is a real challenge to IT professionals to construct a flexible and tractable 
technology to handle large samples of correlated data which will preserve all global properties 
and interconnections of the principal components of the stricture in all data transformations 
and exchanges. 
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